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Measurement of magnetic oscillations on thermodynamic quantities (like magnetization and spe¬ 
cific heat), is one of the experimental methods to access the density of states of electronic systems. 

In the present paper we therefore theoretically explore the oscillatory phenomena on the specific 
heat of graphenes considering gapped and gapless cases in a quantized magnetic field. Further sit¬ 
uations is also considered, as the influence of impurities. Coulomb interaction and phonons. We 
could then map the magnetic oscillations on the specific heat of graphenes under these constraints 
and the obtained results are a good starting point and guide for further experimental works. 
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The band structure of monolayer graphenes shows that 
it is a gapless semiconductor with a linear energy spec¬ 
trum near the so-called Dirac points^ and, as one of 
many consequences, Landau levels (LLs) are non equidis¬ 
tant; and therefore unusual behaviors upon a magnetic 
field change is observed^. In addition, the distance be¬ 
tween Landau levels is very large (for example, the energy 
gap between the first two LLs in a magnetic field of 10 
T is more than 1000 K). This leads to giant magneto- 
optical^i^ and thermo-magnetic effects^"— , as well as to 
an unusual quantum Hall effect, which can be observed 
even at room temperatureiii^”— . 

Due to the valence of carbon ions and the gapless en¬ 
ergy spectrum, the electronic density of states (DOS) of 
graphenes vanishes at the Fermi energy and therefore any 
further knowledge on this issue is useful for the scientific 
community. In this direction, experimental access to the 
DOS could be obtained by measurements of magnetic 
oscillations of the heat capacity under a quantizing mag¬ 
netic fieldi^i^. 

Thus, the present effort theoretically explores the 
oscillatory behavior of the specific heat of monolayer 
graphenes under a constant area, considering the gap¬ 
less and gapped cases; and, in addition, the influence of 
impurities. Coulomb interaction and phonons. The ob¬ 
tained results are therefore a good start point for further 
experimental studies. 


II. MAGNETO-OSCILLATIONS OF 
ELECTRONIC SPECIFIC HEAT 


The thermodynamic potential of the two-dimensional 
electron gas in a quantized magnetic field can written 


Q{T) = - 


2eHkBT§ 

IT he 


El- 


1 + exp 


ksT ) 


(1) 


where H is the magnetic field, e represents the electron 
charge, c = 3 x 10^^ cm/s the light speed, S stands for 
the graphene area and // > 0 represents the chemical po¬ 
tential. From now on, we will use the quasi-classical ap¬ 
proach, based on the quantization conditions of Lifshitz- 
Onsager— . Considering the plane of the electronic sys¬ 
tem along xy directions and an ortogonal magnetic field 
(along 2 : axis), this approach starts considering a quan¬ 
tized area enclosed by an electron trajectory in momen¬ 
tum space: 


A(e) 


2'KheH . 

-(j +7±) 


( 2 ) 


Above, 


777 * 

7 ±= 7 ±^. ( 3 ) 

where (i) 7 is a constant that defines the energy of the ze¬ 
roth Landau level (it assumes 1/2 for non-relativistic gas 
and zero for graphenes), and (ii) m*/2m is the Zeeman 
splitting. In addition: 


1 dA{e) 


2tt de 


(4) 


is the electron cyclotron mass and m the electron mass. 
Finally, j is an integer and e{p) = e. 

For instance, considering free electrons on graphenes 
we would have A{e) = 7re ‘^and then m* = e/vp. For 
graphenes, the values of chemical potential 0 < ImI < 1 
eV are of most interest, since the energy spectrum is lin¬ 
ear for these values^; and thus, from this interval, m* 
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ranges from zero up to 0.1m. These intervals have an in¬ 
teresting consequence: the absolute value of the Zeeman 
splitting then ranges from zero up to 0.05, and therefore, 
since for the present case we are considering |/i| < 1 eV, 
we can thus neglect the Zeeman splitting. 

It is well known that the thermodynamic potential of 
electrons in a quantizing magnetic field consists of two 
components^i^: a zero-field contribution and a 

field dependent term Thus, the total potential 

reads as: 



11(T) = Oo(T) + 

(5) 

whereS^ 

s 

(6) 

and^ 

_ ‘2m*u)ekBTS ^ 1 cos(A:n7r) 

^ nh ^ k sinh(x/c) 

k=l ^ ' 

(7) 

Above, f{e — li) represents the Fermi-Dirac distribution 
function. 
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, 27r2fcBT 

Xk = k 
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and 

eH 

UJr = - 

(10) 


where this last means the cyclotron frequency. 

From the above, the entropy of the system can be easily 
accessed from S = — and reads therefore as: 


where Equation [TT] can then be rewritten 

as: 


So{T) ^ kB 


27rmpSkBT 

3^2 


(15) 


However, in a quantized magnetic field, the chemical 
potential is an oscillating function of the magnetic field; 
but the total number of electrons N must remains a con¬ 
stant and therefore we can consider: 


SAjcF) 


/ h,t 


(16) 


where the first term was obtained considering the cf 
kBT limit on equation [TTl Considering also this limit for 
the second term we then obtain (see reference [H) : 


A{eF)c 

heH 


= a — arctan 


sin(2a) 
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(17) 


where 


a = 


7r‘^ he 


K 


(18) 


and K = A/'/S 10^^ cm“^ is the electronic concentra¬ 

tion. Note the left side of equation [T7I stands for n7r, on 
equation [51 For gapless monolayer graphene it reads as^ 


A{^f) 



(19) 


while for a gapped monolayer graphene^»^ 


A{eF) = ^{el-A^) ( 20 ) 

Vf 


and, finally, for a 2D non-relativistic electrons gas it reads 
as^ 


A{eF) = 27TmeF ( 21 ) 


So{T) = kE 
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cosh" 
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and 


sr(T) = E ( 12 ) 


^2 sinhx/c 
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where 


L{x) = cothx — 


(13) 


Finally, the specific heat (at constant graphene area S) 
can be determined from the entropy above, considering 
Ca = T{dS/dT)a. Thus, the total area specific heat read 
as 




h? 

where 

T(x/c) = 2 cothx/c — Xk 


k=l 


sinhx/c 


kB ( 22 ) 
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sinh^ Xk 


-h cotm Xk 


(23) 


is the Langevin equation. If ei? ^ kBT then integrand 
in equation [TT] is appreciably different from zero only in 
the vicinity of e = /i ~ e^; and thus we can expand the 
function A{e) near e = cf- This assumption leads to: 

(14) 


It should be noted that we consider only the intraband 
contribution to the specific heat, i.e., we considered only 
temperatures kBT hujc- 

Since there are a cosine function with argument de¬ 
pending on the inverse magnetic field iT, it is expected 
that the specific heat oscillates when represented as a 


A{e) « A{eF) + 27rmJ^(e - e^) 
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function of 1/i^; and this effect is caused by quantum 
oscillations of the density of states. See figure [1] for a 
further understanding, where two cases are depicted: (a) 
gapless and (b) gapped graphenes. This oscillations can 
be seen as follows: magnetic field values in which a cer¬ 
tain LL is almost fully filled (and the immediately above 
LL is still completely empty) promotes a sharp change 
on the density of states and, consequently, on the total 
energy and specific heat. A small increasing of the mag¬ 
netic field thus fills the LL and then the specific heat 
reaches to a minimum - and this fact can be physically 
understood, since completed LLs does not exchange en¬ 
ergy (inter-level transition is not allowed, due to the low 
temperature under consideration; and intra-level tran¬ 
sition is also not allowed, since the level is completely 
filled). 

It is also important to stress a changing on the oscilla¬ 
tion period: 


/ 1 \ ^ 27The 


(24) 


comparing the gapped and gapless cases. See equations 
M and [20l as well as figure [T1 


III. LANDAU LEVEL BROADENING 


Real samples have defects and impurities and there¬ 
fore the scattering by these leads to a broadening of the 
Landau levels (known as Dingle broadening). The ther¬ 
modynamic potential considering this broadening of LLs 
can be written as: 


Q,{T) = -ksT j p(e) In 1 + exp 


de (25) 


In the absence of broadening of the Landau levels, the 
density of states is written as: 

n 

leading therefore to equation [TJ however, to consider the 
the broadening effect, the density of states reads as: 


where 


X 1 r 

“ TT e2 + r2 

(28) 

is the Lorentzian approach and 



(29) 


represents the Gaussian one (the present work considers 
only the Lorentzian approach). Above, T is the scat¬ 
tering energy. A detailed study of magnetic oscillations 




FIG. 1. Oscillatory area specific heat for (a) gapless and (b) 
gapped graphenes. 


in graphenes in the presence of impurities, considering a 
Lorentzian distribution, is detailed discussed in reference 
[^ : and, from this work, the thermodynamic potential 
could be obtained: 




2m*uJckBTS 
irh 




exp 


-2'Kk 


cos(/cn7r) 
fujJc\ sinh(x/c) 
(30) 


Considering ep ^ T, the zero field contribution 
^o,r(^) can be neglected; and, on the other hand, T ^ 
completely damps the oscillations of the field de¬ 
pendent contribution (above equation). Thus, the total 
area specific heat, taking the broadening of the LLs and 
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FIG. 2. Influence of impurities on the oscillatory area specific 
heat for (a) gapless and (b) gapped graphenes. These values 
of r are those that satisfy F <C hujc. 


F ^ hujc into account, reads as: 


Ca,r(T) = 


Tirri'kBT 
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oo 
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exp 


(31) 
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where S^{xk) is given by equation [23l 

Indeed, increasing the scattering parameter F, the os¬ 
cillations are suppressed, as can be clearly seen on figure 
[2](a), for gapless graphenes, and[2fb), for gapped case. 
These results are of easy understand, since impurities 
make broader the LLs, decreasing the density of states. 


IV. INFLUENCE OF COULOMB 
INTERACTION 


Coulomb interaction is a very important problem in 
graphene science, since this interaction is not screened 
due to the two-dimensional character of the system. 
To consider this interactions into the quasi-classical ap¬ 
proach, we need to change the conditions of quantization. 
Strictly speaking, nowadays, the problem of the influence 
of the Coulomb interaction on the magnetic oscillations 
is not solved; however, in some studies^^ — this interac¬ 
tion has been introduced into the function A{e) without 
changing the quantization condition. In this paper, we 
will also use this simple approach. 

Thus, the electronic spectrum of graphene renormal¬ 
ized by the Coulomb interaction reads 


e{p) = UbVFP 




(32) 


where g = is the dielectric constant (sub¬ 

strate plu s graphene - for Si02 it is close to 3.9, see ref¬ 
erence l30|), po = b X 10“^^ erg.s/crr^ and assumes 
for the conduction band and for the valence band. 
It is possible to show^, that the electron trajectory in 
momentum space, considering the first order expansion 
in order to small values of is a circle of radius 


a(e) 


e 

Vf 


I-gin 




and therefore 


^(^^) = ^ 


1-gin 


( 77 )]^ 


(33) 


(34) 


The above quantity, for a gapless graphene, can then 
be placed into equation [8] and the specific heat of equation 
[221 revisited. This result is depicted on figure |3l where 
the expected oscillations in the presence and absence of 
the electron-electron interaction can be compared. Note 
the Coulomb interaction leads to a decreasing of the spe¬ 
cific heat, as well as to an increasing and shifting of the 
oscillations. Apparently, this is due to the fact that the 
Coulomb interaction favors a decreasing of the chemical 
potential^^. 


V. PHONONS CONTRIBUTION 


The total specific heat of the system, in addition to 
the electronic term, also contains a phonon contribution; 
that, on its turn, is usually more pronounced than the 
electronic one. As done above, to access this contribu¬ 
tion, we need the thermodynamic potential, that, for this 
case, reads as^: 


^ph{T) = 


SksT 



uj In 


1 — exp 



duj (35) 
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FIG. 3. Influence of Coulomb interaction on the oscillatory 
area specific heat for gapless graphenes. 


three orders of magnitude bigger). Quantitatively, the 
ratio of the phonon and electronic parts of the specific 
heat agrees with results of reference [32. 



FIG. 4. Electronic plus phonon area specific heat for gapless 
graphenes. Note this case is three orders of magnitude bigger 
than the other cases. 


where uu is the phonon frequency and i^ph is the phonon 
velocity (in graphene i^ph ~ 10^ cm/s)^. Introducing a 
new variable y = huj/ksT and integrating by parts, we 
obtain: 

Q 7^3 poo 2 

that, after a simple integration reads as: 

qc?l3 7^3 

= (3T) 

where (^{z) is the Riemann zeta function. Then, the 
phonon entropy is a straightforward result: 

Sph{T) = (38) 


VI. CONCLUSIONS 

The present work evaluates the electronic specific heat 
under constant area for gapped and gapless graphenes. 
Due to the crossing of the Landau level with the Fermi 
energy of the system, oscillations are found, in analogy 
to the de Haas-van Alphen effect, found on the mag¬ 
netization of these materials^, and also on the magne¬ 
tocaloric potentials^ —. Impurities, Coulomb interaction 
and phonons are also considered and deeply evaluated. 
Considering these oscillations are useful for an experi¬ 
mental access to the Fermi surface of materials^, the 
present results are therefore a good start point to use 
specific heat as a probe of the Fermi surface of the stud¬ 
ied cases. 


and, consequently, the total area specific heat: 
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(39) 
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This final result can then be seen in figure jH Note the 
phonon part of the specific heat is much greater than the 
electronic one (oscillations remain, but the total value is 
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